a^ 



Nonperturbative QCD vacuum and 



g>: colour superconductivity. 



H . 

'^ ■ N.O.Agasian* B.O.Kerbikovj V.I.ShevchenW 

Institute for Theoretical and Experimental Physics 

117218, B.Cheremushkinskaya 25, Moscow, Russia 
(N 
> 

in 

en 
m 
(N 

O 

a\ _ 

^ ! We discuss the possibility of existence of colour superconducting state in 

p ■ real QCD vacuum with nonzero {agGG). We argue, that nonperturbative 

p i' gluonic fields might play a crucial role in colour superconductivity scenario. 

> 

X 



Abstract 



*e-niail: agasyan@vxitep.itep.ru 
^e-mail: borisk@heron.itep.ru 
■l-e-niail: shevchen@heron.itep.ru 



1 Introduction 

The behaviour of QCD at high density has become recently a compelhng 
subject due to (re)discovery of colour superconductivity [||, 0. The essence 
of the phenomenon is the formation of the BCS-type diquark condensate 
at densities exceeding by a factor (2 ^ 3) the normal nuclear density 0. 
Colour superconductivity has been studied within different versions of the 
Nambu-Jona-Lasinio - type model |l|, ^, || or the instanton model [0, §. 
To our knowledge, however, no attention has been paid to the fact, that 
as compared to real QCD, both approaches miss important nonperturbative 
gluonic content of the theory. The NJL model contains the gluon degrees 
of freedom in a very implicit way: it is argued, that high-frequency mode 
(one-gluon exchange) after being integrated out gives rise to the effective 
four quark interaction (see, e.g. 0). The instanton model deals only with 
specific field configurations - instantons and antiinstantons. The NJL model 
with gluon condensate included has been considered in [^ while confining 
background superimposed on instantons has been treated in ^j 

The role played by the gluon degrees of freedom in the problem under 
consideration is essentially twofold. First, it is assumed, that they are re- 
sponsible for producing the quark-quark attractive interaction leading at 
small enough temperatures to the Cooper pairing (thus gluons play the role 
of phonons over ion lattice speaking in condensed matter terms). On the 
other hand, the vacuum gluon field fluctuations should be affected by colour 
superconducting state itself in a way analogous to the Meissner effect in ordi- 
nary superconductor. The crucial point is which force will win, i.e. whether 
superconductivity will survive or will be destroyed by gluonic fields as it hap- 
pens in the standard BCS theory in the presence of strong enough external 
magnetic field. 

The studies of idealized QCD performed by several authors have shown. 



that depending on the chemical potential fi and the temperature T the sys- 
tem displays three possible phases: 1) chiral symmetry breaking without 
diquark condensation; 2) mixed phase with nonzero values of both chiral and 
diquark condensates; 3) diquark condensation without symmetry breaking. 
The system may be described by the thermodynamic potential fi(0, A; /i, T), 
where and A are order parameters, related to the chiral and diquark con- 
densates respectively (explicit definitions will be given below). The potential 
Q is expressed in terms of quantum effective action F as fi = TT/Vs. Sup- 
pose, that the potential Q has been calculated within the framework of some 
NJL-type model, i.e. with the gluon sector excluded (except for contribu- 
tion of high-frequency modes, giving the necessary attraction). Consider 
for simplicity the phase 3) of the system with = 0, A 7^ and T = 0. 
Let Q{(f) = 0, Aq; /U, 0) be the stationary value of the thermodynamic poten- 
tial, where Aq is the solution of the gap equation dQ/dA = (see below). 
Now we superimpose the nonperturative vacuum gluon fields on the above 
picture. The detailed knowledge of the nonabelian Meissner effect is unfor- 
tunately absent. Anyway, it is obvious that the corresponding microscopic 
picture is far from being trivial. The nonlinear character of the equations 
of motion for gluon fields is of prime importance here, while usual Meissner 
effect for abelian fields is essentially linear phenomenon. On the contrary, 
general symmetry arguments tell us, that part of gluon degrees of freedom 
becomes massive if the colour gauge invariance is spontaneously broken. It 
means effective screening of low-frequency modes and therefore it is reason- 
able to assume, that the formation of the colour gauge invariance - breaking 
diquark condensate should lead to the decrease of the gluon condensate by 
some factor, which we assume to be about a few units (but, presumably, not 
to exactly zero value, as it happens in abelian case). Then it will be ener- 
getically favourable for the system to remain in the colour superconducting 



state with A 7^ only if the quantity 
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is less than fi(0 = A = 0; /i, 0) — Q{(f) = 0, Ao;/i, 0). The factor n in (|lD 
represents the unknown rate of decrease of the gluon condensate due to su- 
perconducting diquark state formation. Note, that e(fi; -^ 00) = —e^ac- In 
what follows we will show, that these two energy gaps have the same order 
of magnitude unless the k is close to 1. 

It should be noted, that the analogy with the BCS superconductor in 
external field is somewhat loose here for at least two reasons. First, strong 
nonperturbative gluonic fields are inherent for the QCD vacuum. The consis- 
tent way of analysis should include a set of gap equations for the free energy 
of the system depending on gluon and different types of quark condensates, 
determining energetically best values for all of them simultaneously. The 
second point is the relation between scales, characterizing colour supercon- 
ductivity and nonperturbative gluon fluctuations. In particular, only modes 
with the wave lengths larger than the Cooper pair radius are responsible for 
the supercurrent while the rest do not admit simple interpretation in terms 
of Ginzburg-Landau theory. We leave the analysis of these complicated prob- 
lems for the future. 

Another important remark is in order. It might be naively assumed, that 
if the system under study is in the deconfinement region, the vacuum gluonic 
content may be taken as purely perturbative. There are several reasons, how- 
ever, why it is not the case. The most important one is the following. Finite 
density breaks Euclidean 0(4) rotational invariance and hence chromoelec- 
tric and chromomagnetic components of the correlator {asG"^) ,q enter on 
the different footing. In particular, deconfinement, i.e. zero string tension is 
associated with the vanishing of the electric components, while it is energeti- 
cally favourable for the magnetic ones to stay nonzero (the same phenomenon 



takes place for the temperature phase transition [0). At the same time it is 
just strong magnetic field, which is able to destroy the superconductivity. 

Needless to say, that the colour superconductivity phenomenon is essen- 
tialy nonperturbative effect since it requires the formation of the gap, which, 
being proportional to IVqcd^ cannot be obtained in perturbation theory. 

2 General formalism 

We start with the QCD Euclidean partition function 

Z = fvAV^Vtjj exp{-S) (2) 

where 

S=\J F^^F^J^x + J ^lj{-ti,D^ -im + tfi^,)^ d^x (3) 

We supress colour and flavour indicies and also introduce chemical potential 
fi (only the case Nf = 2 ; N^. = 3 is considered in this paper). Performing 
integration over the gauge flelds one gets effective fermion action in terms of 
cluster expansion 

Z = f VijVijj exp [- f (fx Lo - Seff] (4) 



with Lq = ipi—i'jfidf^ — im + ifi'y4)ip and effective action Sefj = Z) -7 ((^")) 

_ n=2 

where 9 = J d'^x 'ip{x)g'y^A1^{x)t"-ip{x) and double brackets denote irreducible 
cumulants. 

To proceed, one is to make considerable simpliflcation of eq.(^. First, 
only the lowest, four-quark interaction term is usually kept in Sejf- Second, it 
is instructive to consider instead of the original nonlocal kernel some idealized 
local one, respecting the given set of symmetries. In the problem under 
study it is common to choose either instanton-induced four-fermion vertex 
(the choice, adopted in 0, ^) or different versions of the NJL model ^ |, |^ 



including the one, motivated by one gluon exchange |^. In the latter case 
one gets 



S,ff = J d^xd'y {i;{x)^,ty{x)){i;{y)^,ei;{y))D{x - y) (5) 

with D{x — y) = g'^5^^\x — y) and where the coupling constant g'^ with the 
dimension mr"^ is introduced. We assume, as it has already been mentioned, 
that this localized form of the kernel does not mimic all gluonic content of 
the original theory, in other words only some part of gluon degrees of freedom 
participates in the condensate formation. One way to analyse the role, played 
by other ones, would be to consider more realistic nonlocal functions D{x — y) 
(which, in principle, encode all necessary information if we keep only four- 
quark interaction). This will be done elsewhere, while in the present paper 
we work with the local form of the action. Performing colour, flavour and 
Lorentz Fierz transformations and keeping only scalar terms in both ijji/j 
and ipip channels, we arrive at 
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(^(a;)AV'(x))(^(y)A^(l/)) - (^(x)$,^^(i/))(^^(y)$,^(x)) (6) 



where 



A ^ " 1 (2) 

A = -^ Ic Tf , $„ = -j= eo/37 75 t'p 

and ip^ = Cip'^ = 7274 ip'^. We note, that with only scalar terms kept, the 
Lagrangian (|^) is no more chiral invariant. The attraction in scalar colour 
antitriplet channel (which also exists if one starts from the instanton-induced 
interaction) could lead to the formation of the condensate, breaking colour 
SU{3). In close anology with we replace the common coupling constant 
g"^ by two independent constants gf and 5f| corresponding to the two terms 
in (^). Next step is to write down the partition function and to perform its 
bosonization. We adopt the standard Hubbard-Stratonovich trick and get 



Z = fV'ipVip exp{- j (fx (Lo + Leff)) 



/ VcpVAVA^exp I f d'^xl- 
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+ TrLn| . '^^^ ^^ + ^im + ^) - ^,,, 

For the system of the massless quarks at the phase 3) and T = the 
thermodynamic potential reads (we remind, that Z ~ exp{—S) ~ exp{+r)) 

Vt= — = ^-4, Tr Ln . '^ 8 

^4 gl V^ \ id'^ + ^74/i 2AtC-i<l>t ) ^ ' 

The value of the diquark condensate is determined by the gap equation 

= (9) 



dA 



A=Ao 

By going in (j^) to wave number-frequency space and making use of the gap 
equation (|^), we arrive at the following expression for the thermodynamic 
potential at its minimum 

^(0 = A = O;/i,O)-^(0 = O,Ao;^,O)- ^ ^° ^ (10) 

where M is the NJL cutoff which is typically about 0.8 Gev [|T], ^. Alterna- 
tively, the cutoff may enter via the formfactor or the instanton zero-mode. 
In the NJL-type calculations the values of the cutoff M and the coupling 
constant f^^ ^^^ fitted simultaneously, but no unique "standard" fit exists so 
far [^. To estimate the r.h.s. of ( p!OD we have taken for cutoff M = 0.8 Gev, 
for coupling g^ = 12(7| = (15 -ir 40) Gev~'^ and the value of the gap in the 
diquark scalar sector Aq = (0.1 -7-0.15) Gev. With these parameters one gets 

(](0 = A = 0; /i, 0) - fi(0 = 0, Ao; /i, 0) ~ (1 ^ 5) ■ 10"^ Gev^ (11) 

To be on a robust quantitative footing and to consider finite temperatures one 



can replace the estimate (11) by the result of direct numerical calculations 



of the thermodynamic potential performed in [^]. Our result (|TT]) is larger 
than the corresponding value, presented in [Q. The discrepancy may be 
due to different value of the coupling constant g^ adopted in 0. Needless 
to say, that the larger is the estimate of ([10|) the larger is the critical field 
extinguishing superconductivity. 

Now let us estimate the expression (|I|). We use different sets of data from 
||TT| and take for the gluon condensate 
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G2 = {— Gl^Gl^) = (0.014 ^ 0.026) Gev 
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n^ 



Then for two flavours one gets in one loop 

~ ("l - -V (4 - 8) ■ 10-3 Gev^ (12) 

It is seen, that (|10D and ([l^) have the same order of magnitude for k > 2. 

It should be noted, that the estimate (0) given above is rather optimistic 
in the following sense. If one naively assume the colour superconductor to 
be the BCS one and take its typical parameters, for example, from 0, then 
one has at T = for the value of the critical magnetic field 

Hl^, = 0.64 ■ (Ao)2 ■ (m.pp) (13) 

where p'p = fi'^ — rn"^, rrig is the constituent quark mass. We take fi = 
(0.4 -^ 0.5) Gev and corresponding Aq = (0.04 -=- 0.10) Gev from the paper 
PI which are smaller, than the typical values we have analysed before. The 
maximum of ([T^) with respect to pp at the fixed fi is reached for /x = \/2pF 
and it gives 

Hl-^ = (0.8 ^ 8) ■ 10-^ Gev^ (14) 



Strictly speaking we are not allowed to apply formulas like (|T3D to the colour 
fields and interpret them in terms of Meissner effect, but we note, that (p!^ 
even without any numerical factors is about order of magnitude smaller than 
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3 Finite density effects 

The comparison made in the previous section was intuitively based on the 
analogy with the Meissner effect in ordinary superconductors. The actual 
value of the gluon condensate in (|1^) was taken to be the vacuum one. This 
is not quite correct, however. Even without the formation of any diquark 
condensate, the gluon condensate in the hadronic matter is different from 
that in the vacuum. In order to get an idea about such dependence, let us 



consider effective dilaton Lagrangian [12 



L{^) = lid.af - V{a) ; V{a) = ^ a^ (in ^ - i) (15) 

where dilaton field is defined according to 



64|eJ ' ' '^'^^ ' 4a, 



and 



tjiq 2 16|et,| 



A = -—, r , CTn 



16|ej,| uiq 

where e„ = | {6|_l^l) is the nonperturbative vacuum energy density and trq - 
dilaton (i.e. glueball 0"'""'" in our case) mass. Low energy dilaton physics can 
be used for description of the gluon condensate behaviour at finite density 
and temperature |13|. In the chiral limit masses of nucleons in QCD are 



determined by the nucleon-dilaton vertex L^nt = Tn*j^qq with the effective 
mass m*j^ = niiq ■ (ct/ctq). In isotopically symmetric system the energy density 
takes the form: 

e(a„, n) = V{a„) - V{a^) + 2 j -^ sjp^ + m%^ (16) 

where n = Un+Up = 2p|,/37r^ is the baryon density and the chemical potential 
fi^ = Pp + rn]^'^. Being interested in the densities close to the nuclear density 



1/3 

rig =0.1 Gev we minimize the total energy density with respect to the 

dilaton field 

de_ 

and find in the region nrriN <C lQ\ey\ 

{ — G" C )n = { — G1ivGliv)o ■ 1 ~ c — ) (17) 

where the factor c = mArno/4|e^| ^ 0.05. It means, that the decrease of the 
condensate at the densities of interest is about (10 4- 25)%. 

However, we should point out, that the matter may have already the 
structure of cold quark-gluon plasma instead of being hadronic at the den- 
sities of interest. Presumably it would require another type of estimates for 
the value of the gluon condensate. The reader is referred to the paper H 



where it is argued, that the long-wave part of colour magnetic field which can 
lead to the Meissner effect is supressed compared to the critical field by the 
perturbative coupling constant (which is assumed to be small due to asymp- 
totic freedom). It is important to distinguish the quantum fields considered 
in |T^ which propagate in the quark-gluon plasma from the nonperturbative 
fields discussed by us here, hence the magnitude of the latter is an "external 
parameter" with respect to the colour superconductivity problem itself. This 
question needs further investigation. 

4 Discussion 

We have confined ourselves in this letter to a very modest aim - to compare 
the energy gap, typical for the colour superconductor and the contribution to 
the vacuum energy density, coming from the gluon condensate. We argued, 
that if there is a competition between colour antitriplet scalar diquark and 
gluon condensates, which is natural to assume from abelian analogy, then 



10 



it happens on the same energy scale and therefore to take this effect into 
account is important for the selfconsistent picture. It is clear, that further, 
more quantitative analysis is needed. 

The idealized case of two flavours was investigated in the paper. For 
chemical potentials exceeding the strange quark mass the Nf = N^ = 3 
scenario is more physical. It seems, that inclusion of the third flavour does not 
change our analysis crucially. Moreover, due to the phenomenon of colour- 
flavour "locking" d], eight of the nine {{N^ — l)su{3) + ^u(i)) gauge degrees 
of freedom acquires mass due to the Higgs mechanism in Nf = 3 case and 
therefore colour superconductivity is "complete". According to the line of 
reasoning adopted in this paper, it means, that K,{Nf = 3) should be larger 
than ^{Nf = 2) (for the scalar channel). 

To conclude, we have discussed the role played by gluon degrees of free- 
dom in the colour superconductivity and, in particular, we have argued, that 
nonperturbative gluon fluctuations might be strong enough to destroy the 
phenomenon. More quantitative analysis of this problem is in progress now. 
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